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IIWEODyCTION 


The description of the acoustical properties of duct liners 
depends on the measurement of properties in an impedance tube* 
However^ in practical appl ications , these liners are usually 
mounted on a side wall of a flow duct* Thus, both the mean flow 
and acoustic wave qraxe over the surface of the liner. Lester 
and Parrott M,21 find that these acoustical liners often act 
differently in cjra*-*ing incidenco than one would expect from the 
impedance tube measurements . The study of these differences is 
hindered by t!ie lack of a lUmple tochnit|ue to experimentally 
determine the properties of duct liners in grazing flow* 

This paper describes a simple technique to measure the 
acoustical properties of duct liners in qrazinq flow. The tech-- 
nique is developed explicitly for rectanqular ducts though it can 
be applied equally well to circular and annulvnr ducts. The tech-- 
nii|ue is to measure the axial wavenumber of the least attenuated 
mode. This wavenumber depends on the (\c\^ profile in the duct, 
on the bounviary condition, and on ttie acoustical properties of 
the duet liner. Thus, the acoustical properties of the liner can 
be found l^Y moasurinq the wavenumber and knowing the flow profile 
and boundary condition. The techn Ique wi 1 1 be introduced for 
unlfvnm flow ducts with point reacting liners. It is then 
av^pliod to bulk reac'ting liners with both rigiil and limp material 
structures. PinalVy the influence of shear flow is considered. 

Tlve acoustic%-^l properties of materials are usually measured 
in im impedance tube. An Impodauce tube is a hollow duct termin- 
ated by a liuct t i 1 Veil with the unknown acoustical material. A 
sound wave impiniies on t!\e acoustical material and is reflected 
back . by meas ur i ng t he i n to r f erence pa t tern bet va^en t he i nc idon t 
and reflected sound waves the acoustical prov’^ertles of the 
material can be determined. If the material Is locally reacting, 
a single measuroment of the inter ferenv'e pattern is needed to 
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find the admittance at each frequency. If the material is 
nunlocally reacting, two measurements must be made at each 
frequency. Two measurements are needed because bulk reacting 
materials are described by two parameters: the propagation 

•constant and the plane wave impedance. 

Three impedance tube techniques exist for measuring the 
acoustical properties of bulk materials. Scott (3) proposed a 
direct method for cases where the length of the bulk reacting 
material is effectively infinite. It is not precisely an 
impedance tube measurement because the propagation constant is 
determined directly by measuring the phase and amplitude of the 
pressure wave with a microphone traversing the length of the bulk 
material. The plane wave impedance is measured directly by 
simultaneously measuring the acoustical pressure and velocity at 
the surface of the bulk material. Dinardo (4| and Romero (5j 
propose an' impedance tube measurement called the two-cavity 
method. A finite length of bulk material is terminated first by 
a hard wall and then by a quarter wavelength resonator. The 
propagation constant and characteristic impedance of the material 
are deduced from a set of transcendental equations. Ferraro and 
Sacerdole (61 propose an impedance tube measurement known as the 
two-thickness method. Measurements are made using two lengths of 
bulk material terminated by a hard wall. Again, the acoustical 
properties of the bulk reacting material are deduced from a set 
of transcendental equations. 

An experimental comparison of these three techniques is 
given by Smith 171. There are two advantages to the two- 
thickness and two-cavity methods. First they avoid the need for 
measurements in the bulk material. Second, they give the acous- 
tical properties by manipulation of simple transcendental equa- 
tions. Those advantages are retained in the method proposed 
herein. 
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Consider a two-dimensional duct with a hard wall on one side 
and a soft wall treated with an acoustical material on the other 
side. The axial wavenumber for each mode propagating in this 
duct depends on the flow profile, the bound.iry conditions, and 
the acoustical properties of the liner material. Hence, by 
measuring the amplitude and phase of the least attenuated mode 
the liner properties can be deduced. The axial wavenumber can be 
determined frtxn a pressure measurement made on the hard wall side 
of the duct. For point reacting materials a single measurement 
at each frequency Is sufficient. For bulk reacting materials two 
measurements at each frequency are needed. These two measure- 
ments can be done with two thicknesses of the liner material, or 
by placing first a hard wall and then a quarter-wavelength 
resonator on the exterior, or with two duct widths. The 
important aspect is to create conditions so that two different 

values of the axial wavenumber are measured. 

« 

The method will be described in stages. First, It is 
described for point reacting liners used in uniform flow ducts. 
The influence of the boundary condition, the method for the phase 
and attenuation measurement, and special problems associated with 
tlie axial extent of the liner are discussed. Next, it Is applied 
to bulk reacting liners; rigid liners are considered first, limp 
liners are considered next. Finally, the influence of shear in- 
flow is discussed. 

Two descriptions of bulk reacting materials are possible. 

One approach, which we will not use, it to tittempt to describe 
the material roicroscopioally and determine sound propagation from 
first principles. Simplified models such as tubes (81, parallel 
fibers [9,10), and arrays of spheres (111 have been proposed, 
Ingard (12} points out that "If in such studies empirical para- 
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meters are introduced, their value is questionable in comparlRon 
with what can be learned frcsn dimensional analysis or from a 
macroscopic phenomenological approach." 

. This macroscopic approach is described by Ingard [12J, 
Zwicher and Kosten [13], Beranek [14,15], and Tack and Lambert 
[16] . The approach is to write the mass, momentum, and compres- 
sibility equations for the air within the duct liner. These 
equations depend on macroscopic properties such as porosity, 
density, and specific flow-resistivity of the material. A wave 
equation is derived from these fundamental equations. The acous 
tical description of the waves depends on the propagation con- 
stant and characteristic impedance which, in turn, depend on the 
macroscopic properties of the material. The experimental tech- 
nique described herein is one that relates the axial wavenumber 
in the flow duct to the propagation constant and characteristic 
impedance Itj the duct liner. The explicit relation is based on 
the macroscopic description of the liner. However, the method 
can be used equally well if other descriptions of the liner are 
used, provided that they admit wave-like solutions subject to a 
dispersion relation and a characteristic impedance. 
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2, TOE UNIFORM FLCm EWT I5ITH A POUfT REACTING LINER 

Consider a sound wave propagating in a two-diicensional 
duct. Let X be the axial direction and y be the transverse 
direction. At y =* 0 let there be a hard v?all, and at y = 1 let 
there be a point reacting liner characterized by an admittance 
a. Between the wall the motion of the air is described by the 
linearized mass and momentum conservation equations* 

PO §t “ ^ ° 



+ PoV*U = 0 


where 


D 

Dt 


,3_.Mc3 
'3t 


3x‘ 


3p 


= c2 


( 2 ) 

(3) 

(4) 


Here, p p and u are the first order perturberations in pressure 
density and velocity. Me is the mean fJow velocity in the axial 
direction, pq the mean air density, and c is the adiabatic sound 

speed. We will search for propagating waves in the axial direc- 

i(k x-ut) 

tion. Hence, assuming all quantities vary as e gives 


ikpoe{l-K^M)u = VP (5) 


i kd-K f1)p = PocV.u (6) 

A 


which, in turn, gives the wave equation 


+ I (i-K M)2 - Kj 1 p = 0 (7) 

k2 3y2 I ’‘I 
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where Kx * kj^/k. 
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The sol on can be constructed from elementary solutions 

i(k x+k y-wot) 

proportional to e " . Thusp the dlspersioh relation 

is 


(1-K H)2 - (K^^+K 2) = 0 (8) 

X X y 

In particular, we must choose a solution that satisfies the zero 
normal velocity requiroment at y = 0. 


p = 

Po cosk^ye 

X 

(9a> 



ik„x 



u^ cosk.v 

e 

(9 b) 

X 

ox y* 

ik X 


U “ 

. sink y 

X 

e 

(9c) 

y 

oy y^ 




The relation between pg, and is fixed by the momentum 

equation. Thus 


Po 


PGU 


ox 


(l-K^M)Kx 


(10 a) 


Pg 


pcu 


oy 


-i(l-K M)K 
X y 


(10b) 


The boundary condition at the soft wall is the particle 
displacement condition. This condition assumes an infinitesi- 
mally thin layer of still air just above the soft wall. The 
particle displacement (t) across the shear layer separating the 
still and moving air is the same. Here, the subscript f indi- 
cates a quantity in the moving fluid and the subscript w indi- 
cates a quantity in the still air just above the soft wall. 
Thus, 
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Prom the momentum equation it follows that in the moving fluid 

Uyg = -ik(l-K^M)5£ and in the still fluid just above the soft 

wall u = ik5, . Thus, at the wall, (1-K M)u = u In turn, 

yw w ' X yw yf 

the momentum condition gives 





Hence, the admittance condition gives 
pqcu 

a = — I y = 1 
= .. -i , ... 3p/9y I 

(l-Kj^M)^ p 'y = 1 


or 


K tan(K kL) = iCl-K H)^a 

y y ^ 


(12) 


Recall that the dispersion relation gives • Thus^ if 

we could measure and knowing k, L, and M, we could uniquely 
determine the admittance a at any frequency a> = ck. For point 
reacting liners only one measurement needs to be done because 
there is only one unknown, the admittance a, in the character- 
istic equation. 


Now consider the problem of measuring k^^. For the moment, 
assume the duct is infinitely long and only the least attenuated 
wave is propagating. Furthermore, assume that a pressure mea- 
surement can be made at any location x on the hard wall y == 0. 
This is an advantageous location because it is nonobstructive and 
also at a pressure antinodo for all modes. Let the constant 
Po = |p| e^^® in Eq. 9. There are two methods to determine the 
real, k^, and imaginary, k^, portions of the axial wavenumber 
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The first method is to define Pg and Pj^ where 

2 rT 

Pg = ^ P(t) COS utdt = IpI e CQs(kj^x4-(j>0 ) 


Pfa = P(t) sin utdt = |p| e ^ sin(k^x+^o) 


-k . X 


U3) 


It follows that the amplitude and phase are 
A(x) = In (Pg^+Pb^) 

= ln(|p|/p^^f)-k.x 


(14a) 


4>(x) = tan tPb/Pa^ 


(14 b) 


Thus the amplitude and phase vary linearly with distance. Hence, 
kj, and kj^ may be found from the slopes of the experimentaliy 
determined amplitude and phase measurement. 


The second method is to determine the temporal cross- 
correlation between two points separated axially. Let 

12(t) =1 /o Pi (t)p2 (tT )dt 


R] 


(15) 


= ipr 


-k j^(x2-xi ) 


cos(kj,(x2-Xi ) - u)t ) 


Now let Rj 2 be normalized to the self correlation at point 1 to 
find 


— Rl2(t) -kj^(X2“Xl) 

Rl2(t) = • r" ( 'c) J ~ ® cos(kj.(x2-xi )-u)T ) (16) 


Thus kj, is found by noting the phase shift of R |2 by 

noting the amplitude of R 12 . 

Finally, consider the practical problem that the length of 
lined duct is usually finite. There are two problems: the 
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influence of high order modes and the influence of reflections. 
First consider the problem of high order modes. Let the duct be 
made of two parts. On one side is a hard wall duct in which a 
sound source is located. On the other side let one wall of the 
duct be lined. At the juncture there will be a reflected wave 
and a transmitted wave in the lined section (for the moment we 
neglect any reflected waves in the lined section). This trans- 
mitted wave will be composed of higher order modes as well as the 
least attenuated mode. These high order modes decay very 
rapidly. Thus, a plot of amplitude vs distance from the juncture 
will reveal a region of rapid decay, where the high order modes 
are influential, and a limit region controlled by the propagation 
of the least attenuated wave. The least attenuated wavenumber 
can be only be measured accurately in this limit region. How- 
ever, the region is determined experimentally as indicated above. 
Furthermore, the region of influence of the high order modes can 
be diminished by making the duct v/idth smaller. This observation 
follows by noticing in a hardwall duct the high order modes 
attenuate at a rate = i (rut/D^-k^ where L is the duct 

width. In general, only use the technique below the first cutoff 
frequency for a hard wall duct of identical width to the test 
section. 

Now consider the problem of reflectiors in the lined sec- 
tion. Let the duct be made of three parts. On either end is a 
semi-infinite length of a hard wall duet. In the center is a 
limited region of duct having one wall lined. Near each juncture 
there will be a limited region where high order modes are influ- 
ential. This region is determined as before. In addition there 
will be a reflected wave from the second junction. This reflec- 
ted wave will be smaller than the incident wave in the lined sec- 
tion. In addition, it will be attenuated as it travels back to 
the inlet junction. If the liner length is sufficiently long, 
the reflected wave will have no influence on generation of sound 
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in the lined section. In effectf the lined section is semi-- 
inrifixte. To determine experimentally if this condition is 
reachedr one must test 1 o liners of different axial lengths. If 
the measured conditions at the inlet juncture are the same for 
both lengths, then one can conclude that the reflected wave is 
not influential. At some d..stance away from the juncture the two 
measurements will begin to differ. Where they begin to differ is 
the limit of the range in which can be measured accurately. 
This will be some short distance away from the exit juncture of 
the shorter lined duct. This region of influence will be the 
same for both the short and long lined duct. Hence the measure- 
ments using the longer lined duct can be used to determine k^. 

Of ^course it is always possible to solve the equation of 
motion exactly and fit the observations directly to some mathema- 
tical model. However, that approach lacks the simplicity of a 
direct mGasurement of the axial wavenumber in a lined duct con- 
taining only the least attenuated wave. This least attenuated 
wave depends directly on the flow profile in the duct, the boun- 
dary condition at the wall, and the acoustical properties of the 
lined material. Presumably, if a clean measurement is made, the 
experimenter will be better able to determine which of these 
three assumptions is incorrect, should discrepancies arise 
between the predicted and observed values of the axial wave- 
number. 
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3. E»."ERSIOH OF THE MKTHOD TO RIGID PULK REACTING I,IdERS 


Consider the liner by describing its macroscopic prGperti<-s. 
Let the liner be isctropic and let the material struct jre of the 
liner be rigid. Furthermore, assume there is no mean flo^^ in the 
liner. In the lined material, the liiiearized mass and momentum 
conservation equations are 


P^gh 


dll 

3 


4- r 


yp == 


(17) 


i£ 

at 


s + 



0 


(18) 


where 



(191 


Here pg and and are the first order perturbations cf th 3 
pressure density and velocity in the lined material, is tho 
mecn air density r g is the structure factor, b is the porosity, 
and Cg the complex phase velocity. Note that we use the avorag<» 
particle velocity in the porous space and not the volume velo- 

city which is hUg» 'The subscript s will be used to identify 
acoustic qualities in the lined material. The inertial mass den- 
sity is larger than the average mass density p^h by the strv:clu"-e 
factor g which is normally between 1.2 and 2. The porosity h io 
the fraction of the liner not occupied by the material structure 
and is normally between 0.85 and 0.98 for practical acouotical 
materials. The flow resistance r is def .led by the steady state 
relation ru^ - complex phase velocity Cg accounts for 

the nonadiabatie compressiGn of air in the porous material due to 
heat tranfer from the entrapped air to the material structure. 


Me 

Hencer 


will search for propagating waves in 

i ( k 

assuming ail quantities vary as e 


the 

X 


axial direction. 
-a>t) 

gives 
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K 2 

-ikp^cgh -p. Ug » V pg 


Ik gh ~ •• 

PqC s 




( 20 ) 

( 21 ) 


where the nondimens ional quantities K„ and u are given by 


- Cg 2 

= gh (t:^) 


'I = ^ 


( 22 ) 

(23) 


These equations give the wave equation 


^s (K| - K| ) p i 0 
k 2 3 y 2 s sx s 


(24) 


where Kgj^ = kg^^/k. The solution can be constructed from 


elementary solutions proportional to e 
Thus the dispersion relation is 


^ ^ ^sy ^ ~ 


- f'sx" * W' ■ “ 


(25) 


where Kgy = kgy/k. 

The possible values of Kgy are determined by the boundary 
conditions. V7e will show how to deduce the propagation constant 
ko and the plane wave impedance p /(hp eu ) = P^/g from the axial 
wavenumber measured on the hard wall side of the duct. 

There are two boundary conditions. One occurs at the junc- 
ture of the air within the duct and the liner material, the other 
depends on how the liner is terminated in the transverse direc- 
tion. Three terminations are explored in detail. 

The first termination is none at all. That 's to say the 
transverse extent of the liner is so great that the sound only 
propagates away from the juncture with the duct. Here 
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Ps * Pso ® 


I (kgjj X +kc y) 


u = u e 
sx sxo 




u = u « e 
sy syo 


^'’^sx * ^*^syy> 


(26a) 


(26 b) 


(26c) 


The relation between the P , U and U is given by the 

so sxo syo 

momentum equation. Thus 


P so 

p cu gh s sx 

o sxo ^ 


(27 a) 


*^so _ j, j, 

Po'=‘"syo ~ ^ ® 


(271.) 


Now consider a liner of depth A te.rmjnated by a hard wall. 

At the wall y - L + A the normal velocity u must vanish. Thus 

sy 






ik_..x 


“sx = “sxo“°®(’'sy<y'<^ 


u ~ u 
sy syo 


sin(kg (y-(L + A))]e 




(28a) 


(28b) 


(28c) 


where p /pc u is the same as before and 
^so sxo 


p cu 
o syo 


= ““i ^ K^K 


W *"s^sy 


Pinallyr consider a liner of depth A terminated by a quarter 
wavelev\gth resonator. At the termination y ~ L t a the pressure 
must vanish. Thus 


Ps" Pso'"‘" "sy 


‘"sx" 


(30a) 
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tk X 

“sx " "sxo " 

(30b) 

1.1c X 

“sy“ “syo"^°® ® 

(30c) 


where p /pc u is the same as before and 
so sxo 


^so ^ 4 iii 5, j, 
p cu gh s sy 

o syo 


{31) 


Again, the boundary condition at the soft wall is determined 
by the particle displacement condition. Again, we assume an 
infinitesimally thin layer of still air just above the soft wall. 
The particle displacement across this shear layer is as 

given in Eg. 11. Now the displacement in the linear equation is 
given by 



(32) 


thus 

u, ,(x, L-5,t) = h{l-K M) u„„(x,L-6,t) (33) 

y X sy 

To satisfy this requirement we must set Kgj^ 
continuity of pressure requires that 

p{x,L-6,t) = p(x,L+6,t) 

Taking the ratio of Eqs. 33 and 34 gives 

K tan (K kL) = i(l-K M)^ g (35) 

y y ^ 


= K„. In addition, 
(34) 
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and 


value 


q K ;; 


termination 


i^one 


I !! — tan (K k^) 

‘1K% 

-I cot(K^ kA) 

q ‘ ^ 

' a 


hard wall 


(36) 


quarter wnvo- 
lonqth rosonator 


The Mrat two forma of 6 in Kq. 36 are also qlven Tack and 

Lambert 116). Recall that the dlareraon relation qlvea 

. ..... V -I IK i’-K 2 ) » (K - K ‘ ) . Hence, 

K ^ >• (1 ^.vv C.X » 

iuat as before, a meaamement of the axial wavemimber reveals tivo 
liner ohar.-.ct er i at i ca. Hee.nsse va' wish to extract two qu.anti- 
ties, K and q^u’, in qeneral two me.vsurements must be made. 

These two measurements must yield two sufficiently ditterent 
values of K^, which we call and 

Given K the left hand side of equation, 3h is determined. 

The riqht l^and side contains the unknown 6. by takinq the ration 
of two succesvsive measurements of ft the value vc /q ran bi r .\n 
celled and a t ranscendental equation in this unknown value 

venuHs* 

If two duct Widtivs L, amt I., are us.ed, and the liner is 
effectively infinite in the transverse direction, then this 

equation is 


(1-K 


X % 


• tan ll~K kl-i 


K ^ -K 
s X 


( 37 ) 


(1-K * )' tan (1-K *)•''* kt._. 


K 


-K ^ 

X , 
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If the duct width is constant and two liner depths 
&i and terminated by a hard wall are usedi then this equation 
is 

{ 1-K^2^ ^ ^ 5 kl. ) 

(K 2 -k . 2)W2tan(kA (K 2-K ,2)\ 
s X* \ s ; 

St — --■■ ■ 

(K 2 -K 2)1 /2tan/k&(K 2 _K ,2\ 

S \ S x2 / 


If the duct width and liner depth are constantj and the 
liner is terminated first by a hard wall and then by a quarter 
wavelonqth resonator, then this equation is 


(1 


(1 




2 ) l /2 


2 ) 1 /Z 


t.,„P - x?"''') 


<V ' " x /'''" >:ot(kA(K^= - 


(.19) 


In each case the left hand side is known and the riqht hand 

side depends only on the unknown K . Given K , one can back 

s s 

substitute into the oriqinal equation to find ii^/q. 

The most advantaqeous technique is to use two liners that 
are effectively infinite in the transverse direction. The reason 
is that it results in an equation for K that can bo solved alqe~ 
braically. To determine whether or not a liner is effectively 
infinite, simple increase the liner depth until does not 
ehanqe anymore. 
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4. BKTENSIOM OF THE METHOD TO LIMP BULK-EEACTIHG mTEKIALS 

Consider the possibility of motion in the solid material 
which forms the acoustic liner. The bulk mass density of the 
liner p i is 

PI = hpo + (l"h)Pj^ 

(40) 

- (l-h)Pm 

where the subscript n\ refers to the material properties, and the 
subscript s to the properties of the entrapped air. The motion 
of the air and material are coupled through the momentum equation 


3u. 


= Pi 


m 

3t 


r(u -u ) 
' m s 


(41) 



3 U 

PoOh — I + r(u^-u_^) 


(42) 


In the momentui.i equation for the material structure, we have 
explicitly assumed that its motion is due only to resistive 
action of the entrapped air as it moves past the surrounding 
structure at a rate This is not the only possibility. 

The structure also responds to the surrounding pressure gra- 
dients, and to its internal stiffness. Zwicker and Kosten point 
out that even for zero resistance, coupling between the structure 
and entrapped air is possible by virture of the, inertial coupling 
indicated by the g factor. However, we wish only to illustrate 
how additional material effects alter the calculations described 
herein. In keeping with this aim we limit our attention to 
simple resistive coupling. From Eqs. 41 and 42 v;e find 

r 

r-iwp j 

1 

1 -i ( w/wq ) 



1 7 


( 43 ) 
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where oio * (r/pj). At low frequencies w/w 0 ~ 0 and u_/u “ 

That is to say » the material structure moves in phase with the 
acoustic wave. At high frequencies Uj^^/Ug ~ f/-iup| » i.e. « is 
90® out of phase with u^ and its amplitude is in the ratio of 
laQ/m. This rat’o can be reduced to zero as u is increased 
without limit. Thus, if the material properties are known 
beforehand, a decision can be made as to the influence of the 
material motion. 

Given the ratio u /u , we can derive the modified momentum 

m s 

relation for the air within the porous material 
3u 

-7p = pogh jY" 


where a = -iwp i r/(r-iujp i ) . Thus, wherever r appears in previous 
results, 0 must now be used. The major change is in the 
expression for Kg, which is now 


K 2 = 
s 


(1 + 


I a 


kpocgh 


) 


(45) 


Ingard 112) points out that the effective resistance o can be 
described as a flow resistance r in parallel with the structural 
mass density pj. 

There is an additional change that leeds to be made. This 
change is to include the effects of the notion of the structural 
material on the boundary condition. Again, we assume an infini- 
tesimally thin layer of still air at the s©ft wall. Mass 
continuity requires that 


(C = h(K-iJ 

w ra s m 


(46) 
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Noting that 5 /L * u /u in the structural material gives 
ms m s 

5„-hu«)e^ 

where 

1-h 1 

■g S . ■ — , ,1 . ....I,. 

" (l-iw/u)o) 

The change in the displacement condition alters the character- 
istic equation to 

K tan(K kL) = ia-K^M)2(l+e 16 (47) 

y y A 

The expression for 6 remains the same as given in Eq, 36, Thus, 
even if the motion of the material becomes important, the same 
techniques may be used to discover the acoustic properties of the 
material, ‘ 
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5, APPLICIkTION TO SS5EAR FMM PKOPILES 

Consider the effect of shear flow In the duct. Let the mean 
flow In the axial direction depend on the tansverse coord Inater 
V(y). The linearized equations of mass and momentura Gontinuity 
are 


Dt 


p + Po V 


Oft wT U + U V + 

Dt X y 3y 3y 


„ 2 _ ,, + iE = o 

3t “y ^ 3y ® 


(48) 

(49) 

(50) 


where the adiabatic compressibility 3p/3p = c^ holds and the 

D 


convective derivation 


Dt 


3 3 

^ + V - 5 — now depends on the trans- 

d C d X 


verse coordinate y. Assuming all quantities vary as 
i(k^. x-wt) . 
e *' gives 


i k (1-Kj^H) p = poc V • u (51) 

i k poc (1-Kj^t1) + Uy ly V * i kj^ P (52) 

i k Po c (1-K t1) u = (53) 

X y ay 

where M is a function of y. The wave equation can be derived 
using the commutation relation ^ ~ iy It *^y 

the momentum relation Uy = / (i k pqC (1-K^M)). Thus, we 

recover Hersch and Catton’s (17) results 
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0 = iIe + 

3y2 


2 Kjj 


153 l£ + 1,2 

3y 3y ^ 


- Kj] p (54) 


Because this equation is an ordinary differenttai equation we can 
use the Runge Kutta technique to develop a numerical eigenvalue 
equation^ Mungar and Plumblee {181 use this approach to study 
annular ducts, here we apply their technique to twordimensional 
rectangular ducts. At 

2 3H 

a(y) = 3y (5-,a) 


b(y) = (kL)2 [(1 -KjjM) 2 - Kj] 


(55 b) 


and introduce the column vector ? with P and 3P/3Y as elements; 


fj = P 



(56) 


The second order wave equation can be rewritten as two coupled 
first order equations 



(57) 


or, equivalently 

' 3 I = IF) f (58) 

3y 

where [Fl is the matrix operation defined by Eq. (57). The 
generalized Runga Kutta solution is to divide the duct width into 
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n subsections of length q = L/n. The solution at the j+lth point 
is related to that at the jth point by the rule 

f j+1 = fj + I [nti + 2ni2 > amj + im ] (55) 

where • 

rai = F (y., f.) 

n>2 = F (yj+ I , f. + -J-). 

a 

m3 = P (y^ + I , fj + -y) 
mi, = F (yj + f f f j + q m3 ) 

Each step is a linear transform of the form 

fj+1 == iTj j fj (60) 

which can be combined to yield 

[Tj fo ( 61 ) 

j=0 J 

= [T] fo 

where [xj is the transfer matrix. This matrix is uniquely 
defined by the function M(y) and the constants Kj^ and kL. As an 
eigenvalue problem, Kj^ is the unknown parameter which is adjusted 
to satisfy the boundary conditions. In the application proposed 
herein,. Kjj is the measured axial wavenumber which is used to find 
the boundary conditions. This characteristic equation is found 
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by noting Op/3y)/p = Dq at one wall and B|j at the other. Hence, 
eliminating 3p/3y in Eg. (61) gives 

0 » {(Til + Bo Ti 2 ) - (T21 + Bo T22)} Po . , (62) 

For a hard wall at the point y = 0, we have Bo = 0, For the 
soft wall at y = L = nq, we must again assume an infinitesimal 
layer of still air near the wall. Pressure and particle dis- 
placement are maintained across the layer. Prom the particle 
displacement condition, the acoustic velocity in the transverse 
direction is determined as before. For a point reacting liner 
Op/ay)/p s (1-K^M)2 Pqc(u/p), where, at the wall, u^/p is equal 
to the admittance a. For a bulk liner the velocity in the still 
air at the wall must again be related to the velocity in the bulk 

liner by the rule u , = u _ h(l+e), where h is the porosity 

yw ys 1 . a 

and e the factor accoun'c’ng for the motion of the material. 

Hence, the eigenvalue equation is 

, Tn ( o for point reacting liners 

1 _ «= i(l-Kj^M(L))2 { 

( (l+e)6 for bulk reacting liners 

(63) 

The factor & is given in Eq. (36). Note that for uniform mean 
flow the factor 'rn/kTi 2 is given explicitly as Ky tan (KykL). 

For shear flow (Tii/kTi 2 ) must be determined numerically from the 
measured values kL, Kj^, and M(y). 

If the flow is not two-dimensional, then the above 
mathematical technique cannot be used. In particular, if the 
axial velocity depends on both transverse coordinates V(y,z), 
then the ordinary differential equation in y becomes a non- 
separable partial differential equation in y and z. Twq 
additional correction terms appear in Eq. (54); a^p/az^ and 
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(2K /{1-K M)] [OM/az) Op/az)], An experimental test is to 
determine what extent p and M depend on the other transverse 
coordinate. For the pressure measurementr this can be done on 
the side of the duct with a probe flush-mounted on the hard wall. 

An alternate experimental technique is to perform the 
measurement in an annular duct. The factors a(y) and b(y) are 
changed to functions of the radial coordinate as indicated by 
Mungar and Plumblee [18] . The hard wall at y = 0 is replaced by 
an interior circular cylinder at r = tq. The soft wall becomes 
the exterior circular cylinder. For point reacting liners the 
admittance is defined as before. For bulk reacting liners the 
factor 8 may be redefined by replacing the ».xponential and 
trigonometric functions of the transverse coordinate with the 
appropriate Hankel and Bessel functions of the radial coordinate. 
In addition, the test must be done so that there is no azimuthal 
mean flow on pressure variation. 
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CONCLUSION 


A method to determine the acoustical properties of both 
locally and nonlocally reacting acoustical Materials in grazing 
flow is described. The method is to measure the axial waveniimber 
k/x in a flow duct. Details are presented for a two-dimensional 
duct lined on one side with a bulk reacting liner and on the 
other with a hard wall, though the technique also applies in 
cylindered and annular ducts. Because the axial wavenumber can 
be found from pressure measurements on the side of the duct, it 
should be easy to implement in the laboratory. Sample calcula- 
tions have been illustrated for finite thickness liners back?d by 
a hard wall or by a quarter .wavelength resonator (a soft wall) 
and also for the special case when the liner may be regarded as 
infinite. The effects of mean flow and shear flow have also been 
considered. Furthermore, the influence of the motion of the 
material matrix forming the liner has been considered in the 
special case when liner stiffness can be neglected. 

For point reacting liners the acoustical material is char- 
acterized by its admittance. fJecause one parameter describes the 
material, only one measurement must be done at eac-h frequency. 
Bulk reacting liiers are described by the propagation constant 
Kg, and plane wave impedance 2. Because two parameters must be 
determined, two measurements must be made at each frequency. 

These can be done with ducts of two widths, with two thicxnesses 
of liner material backed by a hard wall, by a single thickness of 
liner material backed by a hard wall and quarter wavelength 
resonator, or to change condition by any permutation of the 
above. 

The important point is to provide sufficiently differenc 
experimental conditions so that two significantly different 
values of the axial wavenumber are measured. 
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